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We have measured the local and global statistics of singularity velocity, v, and have related these 
through the spatial correlation function of v. The distribution of v is a mixture of a mesoscopic 
distribution of global change in the speckle pattern and the distribution for v for Gaussian random 
fields. When v is normalized by the standard deviation of the fractional intensity change, probability 
distributions and correlation function of v approach those for random Gaussian fields. These results 
are directly analogous to the statistics of transmitted intensity normalized by the total transmission 
and provide a unified framework for understanding statistics of speckle evolution and intensity. 
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Coherent waves multiply scattered by a disordered 
medium form a speckled intensity pattern built upon a 
network of phase singularities at nulls of intensity. Cur- 
rent flux vortices circulate about these phasjj singularities 
towards which equiphase lines converge. [2, Q The mo- 
tion of these singularities is a sensitive probe of changes 
in the scattering material and can be used to detect ma- 
terial defects and deformation, [3[ to diagnose cardiac 
fibrillation, 0] and to precisely position small particles. 
[j| Field structures within speckle patterns are generic in 
that these arise in a wide variety of circumstances, Q are 
robust under perturbation and are governed by universal 
statistics. HI The scale of the speckle pattern de- 

pends only upon the angular distribution of the scattered 
field. Perhaps because of the universality of the speckle 
pattern, its structural statistics have only recently been 
utilized to investigate the character of the wave within 
the medium. Studies of mesoscopic fluctuations have 
focused on first and second order intensity statistics in 
space, frequency, [l(| and time, [HI and on their re- 
lation to enhanced fluctuations of total transmission 12 1 
and conductance 13| and to Anderson localization. 14 1 



However, recent studies have found that mesoscopic fluc- 
tuations and localization strongly influence the statistics 
of the structure of static [f| and evolving [TBI speckle pat- 
terns. The statistics of various measures of overall change 
of the transmitted speckle pattern as the frequency is 
tuned were found to be governed by a distribution with 
the same functional form as that for total transmission 
and to similarly depend only upon the variance of the dis- 
tribution itself. Among the measures of speckle change 
are the average displacement or speed of singularities, 
and the standard deviations of fractional intensity change 
or phase change within the speckle pattern. The similar- 
ity in these statistics is surprising since they are oppo- 
sitely correlated with resonances with quasimodes of the 
medium. Whereas the average intensity within a speckle 
pattern for a single incident mode a, I a , which is the to- 
tal transmission, peaks on resonance, the speckle change 
is maximized between resonances. On resonance, the 



speckle pattern is essentially that of the resonant mode 
and so changes slowly with frequency shift. 

In this Letter, we report striking parallels between the 
spatial correlation of local speckle pattern change, mea- 
sured by the generalized speed of phase singularities at 
the sample output with frequency shift, v = dr s /dv, 
and intensity, /. We measure the probability distri- 
bution, P(v), of v normalized by its ensemble average, 
v = v/(v), and the spatial cumulant correlation function, 
C v (Ar) — (Sv(t)v(t + At)), where Sv = v — 1, Ar is the 
distance between the two singularitieslocated at r and 
r + Ar. P(v) broadens and the long-range component 
of C v (Ar), k v , increases in the localization transition. 
k v is nearly equal to the variance of a key measure of 
change of the speckle pattern as a whole, which is the 
standard deviation of fractional intensity change within 
the speckle pattern. This is similar to the near equality 
of the degree of intensity correlation k/ to the variance of 
total transmission. When v is normalized by rj, a quan- 
tity which characterizes the whole speckle change, and I 
is normalized by I a , their corresponding statistics closely 
match those for Gaussian random wave fields. 

We measured the microwave field transmitted through 
samples of alumina spheres contained in a 61-cm-long 
copper tube with the diameter of 7.0 cm. The sample 
is composed of 0.95-cm-diameter alumina spheres with 
refractive index 3.14 embedded in Styrofoam shells of re- 
fractive index 1.04 to produce an alumina volume frac- 
tion of 0.068. [l(| The in- and out-of-phase components 
of the transmitted field polarized along a 4-mm-long and 
0.5-mm-diameter wire antenna are measured with use of 
a vector network analyzer. The spatial distribution of 
the transmitted field over a range of frequencies is ob- 
tained by measuring field spectra at each point on a 1- 
mm-square grid over the output surface of the sample. 
Measurements are made over the frequency ranges 14.7- 
15.7 GHz and 10-10.24 GHz in which waves are diffusive 
and localized, respectively. Frequency steps are chosen to 
be approximately 1/7 of the field correlation frequency. 
Measurements were made in 40 and 71 different configu- 
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rations for diffusive and localized waves, respectively. 

In order to accurately determine the positions of phase 
singularities, the 2D sampling theorem is applied to the 
data to reconstruct the speckle patterns on a 50/j.m x 
50/xra grid. The sampling theorem is also used to in- 
terpolate in the frequency domain, so that spectra with 
120kHz and 250kHz frequency steps are obtained for lo- 
calized and diffusive waves, respectively. This allows us 
to accurately locate phase singularities and to measure 
the magnitude v. Velocity distributions for diffusive and 
localized waves are shown in Fig. 1(a). These results are 
compared to simulations for Gaussian random waves gen- 
erated by the superposition of 300 phased plane waves, 
E(x,y,z) = J2i A* exp[i(k x % + k y y + k z z)\. Each of the 
components of the fc-vector and the amplitude Ai are 
drawn from a Gaussian distribution. We noticed, how- 
ever, the simulation result is not in good agreement with 
the theoretical formula, P(v) = ^^+4)2 1 derived from 
Gaussian random waves in Q. The singularity velocity 
in the x — y plane is tracked as z increases. For diffu- 
sive waves, P(v) is close to the results of the simulations, 
while P(v) for localized wave is noticably broader. Mea- 
surements of C v (Ar) for diffusive and localized waves are 
shown in Fig. 1(b). No corresponding theoretical calcula- 
tion has been done so far. Extremely high value of C„(0) 
is consistent with the fact that the velocity of singulari- 
ties diverges as Ar — > where singularities are created or 
annihilated. C v (Ar) falls rapidly with Ar and reaches 
a constant value denote by k„, which is 0.039 for diffusive 
and 0.648 for localized waves, respectively. Thus the cor- 
relation function can be expressed as the sum of a short- 
range term and a constant: C„(Ar) = C ViS hort(Ar) + k v . 

The statistics of singularity velocity can be compared 
to first and second order statistics of polarized intensity 
which are plotted in Fig. 2(a) and (b). Fluctuations of I 
are seen greatly enhanced for localized waves. The struc- 
ture of C v (Ar) and of the cumulant correlation function 
for intensity, C/(Ar), which is plotted in Fig. 2(b) using 
the same data used in Fig. 1 are similar in that C/(Ar) 
may be expressed as, C/(Ar) = C j tS hort{Ar) + kj. Here 
kj is the degree of correlation, which is the value of 
C/(Ar) at points at which the field correlation functions 
vanishes. [l7[ kj = 0.12 for the diffusive and 3.0 for lo- 
calized waves. In quasi-lD samples with a large number 
of transverse modes, the field in individual speckle pat- 
terns can be assumed to be a Gaussian random variable. 
Thus the probability distribution of intensity normal- 
ized by the average intensity within the speckle pattern, 
I'(r) = I(r)/I a , should be statistically independent of 
the total transmission and the polarized intensity should 
follow the probability distribution P(I') = exp(— /'). 
The measured P(I'), however, deviates from this pre- 
diction [see Fig. 2(c)]. Such deviations are expected 
because the number of transverse waveguide modes is 
small; approximately 30 at 10GHz and 50 at 15 GHz. 
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FIG. 1: (Color Online) (a) The probability distributions and 
(b) cumulant spatial correlation functions of velocity of phase 
singularities normalized to the respective ensemble averages, 
v, for diffusive and localized waves. 



Agreement with Gaussian statistics is better for diffu- 
sive waves since the number of modes is larger. If we 
assume I 1 and I a are statistically indepentent, we find 
(I'(r)I'(r + Ar)/ 2 } = (I'(r)I'(r + Ar))(I%). The cumu- 
lant intensity correlation function, Cp — T/(Ar) — (I) 2 , 
can then be expressed as, 



C/(Ar) = [C,,(Ar) + l] var(/ a ) + 1 



1 



(1) 



Ci' (Ar) and the square of the corresponding field corre- 
lation function, F(Ar), are calculated for diffusive and 
localized waves and seen to be similar in Fig. 2(d). Since 
Cj'(Ar) — > for large Ar, Eq.([T]) gives, m = var(/ a ). 
This is roughly consistent with the measured values of 
var(/ a ) of 0.14 for diffusive and 3.3 for localized waves 
(see Table 1). Again, Cj'(Ar) is closer to _F(Ar) for dif- 
fusive waves rather than for localized waves showing the 
corresponding field statistics is closer to Gaussian when 
the number of transverse modes is larger. 

The character of speckle pattern change can also be 
traced to the combined factors of long-range correlation 
in speckle change and the statistical independence of lo- 
cal and global fluctuations. We consider the motion of 
phase singularities and assume that the global change of 
speckle patterns, denoted by i], is statistically indepen- 
dent of the velocity of individual singularities normalized 
by this change, v' — v/r/. C„(Ar) can then be expressed 
as, 



C„(Ar) = [CV(Ar) + 1] [var(r)) + 1] - 1 



(2) 



If the local changes in the speckle patterns were a Gaus- 
sian random process, we would expect that C v > ( Ar) — > 
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FIG. 2: (Color Online) The first and second order statistics of 
/ and I' for diffusive and localized waves, (a) Probability dis- 
tributions of I. (b) Cumulant correlation functions of I. The 
saparations Ar are normalized by the corresponding correla- 
tion lengths L c , which are the first zeros of the real part of the 
field correlation functions, (c) Probability distributions of /'. 
(d) Cumulant correlation functions of I' and their comparison 
to the square of the field correlation functions. 



for large Ar, since the fields in distant regions would not 
be correlated. In the limit of large Ar, this would give 
k v = v&r(fj). 

Several candidates for a paremeter 77 to quantify the 
change of the speckle pattern were discussed in [15J, in- 
cluding the average velocity of phase singularities v a , 
the standard deviation of phase changes oa^ and the 
standard deviation of fractional intensity change, <xa/» , 
where I* = yfe^arferf^j • ^ ne spectra of these quantities 
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var(7 Q ) 




var(S a ) 


var(ffA<p) 


var(o- A /*) 


var(t£) 


Diff 0.12 


0.14 


0.039 


0.193 


0.087 


0.045 


0.131 


Loc 3.0 


3.3 


0.648 


1.240 


0.743 


0.586 


0.267 



TABLE I: Comparison of ki and k v to variance of global 
measures of intensity and speckle changes. 



in specific sample configurations are similar for localized 
waves and their probability distributions have the same 
functional form as the probability distribution of total 
transmission, though the variances of the distributions 
151 ] The measured variances of v a , cta v and ctaj* 



differ. 

together with n v for both diffusive and localized waves 
are given in Table 1. We find that var(dA/*) is clos- 
est to k v in both cases, while var({i a ) and var(<7A ¥) ) are 
higher. This is also reflected from considerably smoother- 
spectra of CTA7 
been shown in 



than those for v a and cta^, which have 

does 
Thus, 
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In the limit Av -> 0, AI* 
not diverge near singularities as does Aip/ Ai 
var((7Aj* ) more reliably reflects the change of the speckle 
pattern as a whole, while oa v is strongely effected by 
the immediate region around the singularity and v a is 
a property of the small numbers of singularities. This 
suggests that rj = ctaj* is a more practical choice as an 
indicator of global speckle change. 

When v is normalized by 77 = <tai* first and second 
order velocity statistics for diffusive and localized waves 
collapse to the results found in simulations for Gasussian 
waves (Fig. 3), just as was found for statistics of the 
intensity normalized by I a , I'. Despite many similari- 
ties, there are key differences in the relationships between 
global and local statistics for v and /. Unlike intensity, 
which is defined at all points, there are only a small num- 
ber of singularities in the speckle pattern; 15 for diffusive 
and 10 for localized waves on average. This difference 
leads to another dissimilarity between the statistics of 
v and /, that var(/ a ) is slightly larger than kj, whereas, 
var( , D a ) is significantly greater than k v (see Table 1). The 
source of this difference can be seen by first considering 
the relationship of var(/ Q ) and ki. The variance of to- 
tal transmission can be expressed in terms of the spatial 
correlation function of intensity, 



var(/„) - jJ^diA^dAr 2 



Aj A 
= -t 

A J A, short 
~ F J short + K lt 



C^ short (Ar)dAr 2 



+ ki (3) 



where, A is the total area of the output surface. The 
small excess of var(J a ) over k/, indicates that the as- 
sumptions made in Eq. (1) are not strictly valid. A 
quantitative measure of the breakdown of independence 
of II and I a is the relative magnitudes of the contribu- 
tions to var(/ a ) by the integral of the short-range corre- 
lation function C/. s /j 0r t(Ar) over the output surface and 
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kj. Since C/. s hort(Ar) falls rapidly to for Ar > Lc and 
the correlation length Lc is much smaller than the diam- 
eter of the sample cross-section, the integral over A, giv- 
ing Ti^hort = 0.015 for diffusive and 0.088 for localized 
waves, is significantly smaller than the corresponding val- 
ues of he j . Thus the assumptions made are approximately 
valid and var(J a ) « Ki, as expected from Eq.(l). Using 
Eq. (1), we can approximate Tj tS h or t as, 

Ti, short w (1 + «/)- J C r ^ hort {Ar)d 2 {Ar), (4) 

in which, J. C ii tS hort{Ar)d 2 (Ar) corresponds to 
purely Gaussian random fluctuation. 

Equation (3) cannot be applied directly to w because 
the singularities do not exist at every point as does the 
intensity. Finding singularities separated by Ar must be 
described as a correlated random process with a prob- 
ability which is not uniform in Ar. Thus var(u a ) can- 
not be expressed simply as a two-dimensional integral 
of C v (Ar) as was the case for intensity. However, the 
short-range contribution to vax(v a ) can be evaluated us- 
ing the measured values of v&r(v' a ), which corresponds to 
fluctuations for Gaussian waves, multiplied by the meso- 
scopic enhancement factor (1 + k v ) as in Eq. (4) for the 
intensity. We then expect that 

var(£ a ) w (1 + K v )var(v' a ) + k v . (5) 

In the limit k v — > 0, var(£> a ) reduces to the Gaussian 
term v&r(v' a ). Using values of var(u^) in Table 1, Eq. (5) 
gives var(u a ) ss 0.175 and 1.09, for diffusive and local- 
ized waves, which are in reasonable agreement with the 
corresponding measured values, 0.193 and 1.24. 

In conclusion, we have demonstrated a unified frame- 
work for the statistics of transmission and speckle change. 
In each case, mesoscopic fluctuations disappear when the 
local variable is normalized by a global variable reflecting 
the speckle pattern as a whole. In the limit in which a 
large number of modes contribute to the field, first and 
second order statistics of the normalized local variable 
approach the statistics of a Gaussian random process. 
We expect that the statistics of change in the speckle 
pattern, which may arise from internal motion of the 
sample, temperature change, time delay following pulsed 
excitation or by non-monochromatic excitation can be 
described within this framework. 

We thank Bing Hu for contributions to the experiment 
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under Grant No. DMR-0538350. 
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FIG. 3: (Color Online) Probability distributions (a) and the 
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